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Abstract. The main objective of this article is to study the order-disorder 
phase transition and pattern formation for systems with long-range repulsive 
interactions. The main focus is on the Cahn-Hilliard model with a nonlocal 
term in the corresponding energy functional, representing the long-range re- 
pulsive interaction. First, we show that as soon as the linear problem loses 
stability, the system always undergoes a dynamic transition to one of the three 
types, forming different patterns/structures. The types of transition are then 
dictated by a nondimensional parameter, measuring the interactions between 
the long-range repulsive term and the quadratic and cubic nonlincarities in 
the model. The derived explicit form of this parameter offers precise informa- 
tion for the phase diagrams. Second, we obtain a novel and explicit pattern 
selection mechanism associated with the competition between the long-range 
repulsive interaction and the short-range attractive interactions. In particular, 
the hexagonal pattern is unique to the long-range interaction, and is associ- 
ated with a novel two-dimensional reduced transition equations on the center 
manifold generated by the unstable modes, consisting of (degenerate) qua- 
dratic terms and non-degenerate cubic terms. Finally, explicit information on 
the metastability and basin of attraction of different disordered/ordered states 
and patterns are derived as well. 
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1. Introduction 

Many systems in nature can be modeled through the inclusion of long-range in- 
teractions, examples including uniaxial ferromagnetic films, Langmuir monolayers, 
block copolymers, and cholesteric liquid crystals. New features appear in the phase 
transition in such systems where long-range repulsive interactions compete with the 
short-range attractive interactions. Competing interactions can lead to the emer- 
gence of modulated phases, where a particular pattern and wavelength are selected 
[9, 26]. In this paper, we will concentrate on the effects of such an interaction on 
the phase transition in the diblock copolymer melts. 

A diblock copolymer is a linear chain macromolecule composed of two subchains 
of chemically distinct repeat monomers joined covalently to each other. One of 
the subchains is formed of type-A monomers, and the other of type-B monomers. 
Diblock copolymer melts have been studied for decades by polymer chemists due 
in large part to their ability to assemble into various ordered structures at low 
temperatures. However, unlike many other binary systems, e.g. polymer blends 
or polymer solutions, in which macroscopic transition occurs, the phase transition 
here is on a mesoscopic scale where the microdomains of A-rich and B-rich regions 
emerge, and the domains observed have highly regular periodic structures such as 
lamellae (LAM), hexagonally packed cylinders (HPC), spheres, double bicontinuous 
gyroid (DG), and Fddd patterns [2, 3, 14]. 

Through decades of experimental and theoretical investigations, a number of 
elegant theories have been developed to describe the behavior of diblock copolymer 
melts. Most of them are based on mean field approximation; see among many 
others [1, 2, 3, 7, 11, 14, 15, 24]. Also, the classical phase diagram predicted by 
mean field theories or constructed by experiments indicates several regimes for the 
phase separation; see [20, 11, 16] and references therein. 

The main objective of this article is to study the phase transition and pattern 
formation associated with the order-disorder transition (ODT) of noncrystalline 
diblock copolymer melts and to analyze the effects of the long-range repulsive in- 
teractions on the phase transitions. 

The modeling is based on the following considerations. First, for simplicity, we 
neglect the polydispersity effect, i.e., each polymer in the melt has the same poly- 
merization number N and the same number of type-A monomers. Also, we assume 
that both types of monomers share the same Kuhn statistical length I, which mea- 
sures the statistical distance between two adjacent monomers in the same subchain. 
The model we adopt is based on a density functional theory (DFT), proposed by 
Ohta and Kawasaki [23] . The DFT leads to the free energy functional with a nonlo- 
cal term measuring the long-range interaction. Using the time-dependent Ginzburg- 
Landau theory for equilibrium phase transition [19, 17], a modified Cahn-Hilliard 
equation revealing the ODT is obtained, and is also called Ohta-Kawasaki model. 

The mathematical analysis of the model is carried out using the dynamic tran- 
sition theory developed recently by Ma and Wang [17, 18]. The main philosophy 
of this theory is to search for the full set of transition states, giving a complete 
characterization on stability and transition. The set of transition states is often 
represented by a local attractor. Following this philosophy, the dynamic transi- 
tion theory is developed to identify the transition states and to classify them both 
dynamically and physically. One important ingredient of the theory is the introduc- 
tion of a new classification scheme of transitions, with which phase transitions are 
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classified into three types: Type-I, Type-II and Type-Ill. In more mathematically 
intuitive terms, they are called continuous, jump and mixed transitions respectively. 
Basically, as the control parameter passes the critical threshold, the transition states 
stay in a close neighborhood of the basic state for a Type-I transition, are outside of 
a neighborhood of the basic state for a Type-II (jump) transition. For the Type-Ill 
transition, a neighborhood is divided into two open regions with a Type-I transition 
in one region, and a Type-II transition in the other region. 

We now describe briefly the main results obtained in this article. Here we indicate 
only briefly the physical significance of the results, and more detailed physical 
conclusions will be addressed in a separate article. 

First, we show that as soon as the linear problem loses stability, the system 
always undergoes a dynamic transition to one of the three types, forming different 
patterns/structure. The types of transition are then dictated by a nondimensional 
parameter, measuring the interactions between the long-range repulsive term and 
the quadratic and cubic nonlinearities in the model. For example, in the LAM and 
HPC cases, this parameter is given by 

see (4.1) and (4.3) for the exact formula. B > indicates first order transitions, 
and B < implies second-order transitions. 

Second, the long-range interaction (au term in the equation) plays an essential 
role in pattern selection. In particular, there are three unique features caused by 
this long-range interaction term, which are not present in the phase transition dy- 
namics and pattern formations described by the classical Cahn-Hilliard model [19]. 
The first one is a novel pattern selection mechanism associated with the compe- 
tition between the long-range repulsive interaction and the short-range attractive 
interactions, which leads to patterns like LAM, HPC, etc, as mentioned before. 
The second unique feature is that the scale of the spatial patterns emerging from 
the transition have nanoscale rather than macroscale periods, determined by the 
parameter a, which measures the strength of the long-range interaction. Finally, 
the long-range interaction causes the periodic structure of the transition solutions 
to be very sensitive to the system parameters, which can be seen from a precise 
pattern selection criterion (3.6). 

Third, an important technical ingredient of the study is the reduction of the 
partial differential equation model to the center manifold generated by the first 
unstable modes. Then the reduced transition equations are a finite-dimensional 
dynamic system, and are analyzed carefully following the idea and philosophy of the 
dynamic transition theory. Different from many other dynamic transition problems 
we have encountered, in the HPC case, the reduced transition equations are a 
two-dimensional dynamical system, consisting of (degenerate) quadratic terms and 
non-degenerate cubic terms. This unique feature of the reduced system is caused 
directly by the introduction of the long-range interaction term in the model. 

Fourth, in the HPC case, the dynamic transition is a mixed type, Type-Ill. 
The HPC structure is represented by the local attractor for the continuous tran- 
sition part of the Type-Ill transition. Furthermore, these hexagonal patterns are 
metastable. Namely, the original system undergoes a dynamic transition either to 
the hexagonal structure or to some more complicated far away patterns depicted 
by the jump transition part of the Type-Ill transition. 
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This article is organized as follows. In Section 2, a Cahn-Hilliard type equation 
incorporating long range interactions is derived based on the DFT theory. The 
linear problem of this equation is presented in Section 3, and the phase transition 
and pattern formation are studied in Section 4. Section 5 is devoted to the proofs 
of the main theorems. 

2. A Dynamic Phase Transition Model for Diblock Copolymer Melts 

Let ua and be the fraction of monomer number densities of type A and type 
B respectively. Since copolymer melt is highly incompressible, we have ua+ub = 1- 
In the disordered state, ua — a is a constant. The free energy functional in terms 
of Ua derived from the DFT has the following form: 



(2.1) F{ua) = / ^\Wu A \ 2 + f(u A ) 



+ y(-A)"5 {u A - a) ■ (-A)"* (u A -a)dx + F . 



( 2 ' 2 ) ^ ~ ToTTi Z\> CJ d, ~ ]2 V 2„2, , 



Here F is a constant, and the coefficients \x and cr d are: 

kTl 2 _ 9kT 

12a(l - a) ' ° d ~ l 2 N 2 a 2 (l-a) 2 

where I is the Kuhn statistical length as mentioned before, k is the Boltzman 
constant, and N is the total number of monomers in each diblock copolymer chain. 

For the energy term /(ua), we use the following Flory-Huggins formulation for 
the Gibbs free energy of mixing: 1 

(2.3) f{u A ) = kT[^- lnu4 + 1 " A ln(l - u A ) + X«a(1 - u A )}, 

i\a w b 

where Na and Ng are the number of typc-A and type-B monomers in a diblock 
copolymer chain respectively, and \ is the Flory-Huggins parameter which measures 
the incompatibility between type-A and type-B monomers. This expression was 
first derived to describe the free energy of polymer solutions and polymer blends 
[8, 10, 13, 27]. 

The first term in (2.1) represents the interfacial free energy, the second one 
measures the bulk energy of the mixing, and the third term reflects the long range 
interaction caused by the connectivity of the monomers in a chain, where (— A) -1 / 2 
is a fractional power of the Laplace operator under zero flux boundary condition; 
see [6] for more details. 

We note that the original formula derived in Ohta and Kawasaki [23] was given 
for the whole space, the expression on a bounded domain first appeared in Nishiura 
and Ohnishi [21], and the exact expression of the coefficients that we adopt here 
follows Choksi and Ren [6]. 

The equation governing the dynamical behavior of the model can be derived as 
follows. 

Let J be the flux of type-A monomers. Then 

J = -mV(p A - Ms), 



4n [6, 21], the energy term /(«a) (W(ua) in their notation) is a double well function with 

2 a( 



global minima at ua = and ua = 1, typically of the form 4,u a (ua — l) 2 - The variable is 



u = ua — ub- Under this notation, the typical expression of W(u) is l/4(« 2 — l) 2 , which is the 
one given in [21], 
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where m is the mobility measuring the strength of diffusion (see [22, 4]), and /ia 
and hb are the chemical potentials per monomer of Type-A and Type-B monomers 
respectively, which satisfy 

HA- V-B = ~? F(ua), 

where SF(ua)/Sua is the variational derivative of F(ua)- See Reichl [25], Novich- 
Cohen and Segel [22] and Ma and Wang [19]. 
By conservation of mass, we obtain 

We assume for simplicity that m is a positive constant. Then (2.4) is reduced to 

Now let u = ua — a. Using (2.1) and approximating the / term in (2.3) by 
its Taylor expansion about u = up to the fourth order, we obtain the following 
equation for u from (2.5): 

^ = m[-^A 2 u + A(biu + b 2 u 2 + b 3 u 3 ) - a d u], 

(2.6) 



/ u(x, t) Ax = 0, 
Jn 

where bi, b 2 , and b 3 are as follows: 

2 + (l 



(2.7) b 2 = 



a 2 (l-a) 2 
kT a 3 - (1 - a) 3 
2N a 3 (l-a) 3 : 
kT a 4 + (1 - a) 4 
W a 4 (l-a) 4 ■ 



2 X N 



The second equation in (2.6) is due to the fact that the copolymer melt system 
has no material exchange with the external environment. The equation (2.6) is 
supplemented with Neumann and no-flux boundary conditions: 

du dAu „ 0/ _ 
2.8 — = — - on dn, 

on on 

and for simplicity, we consider the spatial domain to be a bounded rectangular 
domain, i.e. fl = lf l=1 (0, Lf). 

We start with the non-dimensional form of (2.6). Let 

d 4 , , , d 2 6i 

x = ax , t = 1 , u = u , A = , 



m/i 

d 2 b 2 d 2 b 3 

) 73 = 1 G d=-ji 

li H a 
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where 
have: 

(2.9) 
(2.10) 
(2.11) 
(2.12) 



d is a typical length scale related to the domain Q. By (2.2) and (2.7), we 
\2d 2 a{\ -a) 



A = 



I 2 



_ 6d 2 [a 3 - (1 - a) 3 } 

72 ~ l 2 Na 2 (l-a) 2 ' 
_ 4rf 2 [a 4 + (l-a) 4 ] 

73 ~ l 2 Na 3 (l-a) 3 ' 

108d 4 



2v a 2 + (i-q) 2 

X Na 2 (l-a) 2 



a = 



l A N 2 a(l-a)' 



With the above non-dimensional variables and parameters, the equation (2.6) 
together with the Neumann and no-flux boundary conditions and initial condition 
can be rewritten as follows (omitting the primes): 

-A 2 u - XAu + A(7 2 u 2 + 73U 3 ) - cm, 



du 

dt 



(2.13) 



/ u(x, t) dx = 0, 
Jn 

du dAu 

— = =0 on dil, 

on on 

u(0) = V, 



where ft = IlLi( > L i)> and L i = L fM 1 < « < 3. 
For the mathematical set-up, let 



, dx = 



H= jue L 2 (ft) ^ 

on on J 



We define the operators L\ = —A + B\ and G : Hi — » by 

Am = A 2 u, 
B a m = — AA?i — cm, 
G{u) = -f 2 Au 2 + 73 Au 3 . 

Thus, equation (2.13) is equivalent to the following operator equation: 



(2.14) 



(2.15) 



du 

— = L A w + G(w), 
it(0) = V- 



It is well known that the Cahn-Hilliard equation can be viewed as a gradient flow 
with respect to H -1 . This is also true for our equation (2.13). Thus the existence 
and uniqueness of solution to (2.13) and existence of global attractors can be proven 
in a standard fashion. 
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3. Principle of Exchange of Stability 

In this section, we consider the eigenvalue problem of the linear part of equation 
(2.15). Let pk, &k be the eigenvalues and eigenfunctions of the following eigenvalue 
problem: 

- Ae K = p K eK, 
de K 

(3.1) ^T = ° onao < 

/ exdx = 0. 

Because of the Neumann boundary condition, px and ex must have the following 
form: 

PK = \K\\ 

(3-2) 

e K = cos( — — ) cos( — — ) cos(— — ), 

L\ L12 

where K is chosen from the following permissible set V: 

(3.3) P=|(^ r ,^ r ,^ r )|fc, e Z-,l< i <3,gfc^o|, 

and 

(3-4) \K\ 2 = n 2 f2kVLl 

i=l 

It is clear that the linear operator L\ = —A + B\ defined by (2.14) with (2.8) 
has the same eigenfunctions {e^ | K e V} as in (3.2). The eigenvalue of L\ 
corresponding to each ejf is given by 

(3.5) (3 K (X) = -|^| 4 + \\K\ 2 o = \K\\\ i^p^)- 

The linear stability and instability are precisely determined by the critical- 
crossing of the first eigenvalue, which is often called principle of exchange of stabil- 
ities (PES). For this purpose, we define a number A c by the following: 

(3-6) ^ = S^' 

Now, let 

(3.7) S = { K e V | K achieves the minimum in (3.6) } , 

then the following PES condition holds true: 
< if A < A c , 

= if A = A c , for K e S, 

> if A > Ac 

(3.9) Pk(Xc) < for all K e V\S. 

The above PES condition shows that A c is exactly the critical value of A, such 
that the disordered state u = is linearly stable when A < A c , and becomes linearly 
unstable when A crosses A r from below. 



(3.8) (3 K (\) { 
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Note that A c is achieved at some K c e S such that \K C \ 2 w y/a. By (2.12), this 
can be written in the dimensional form as: 

1/2 



(3.10) \K C \ = 



^ {Lj) 2 J /7iV</a(l- a ) : 



where £^(1 < i < 3) is the dimensional size of the domain. Since the Kuhn 
statistical length I is of molecular length scale, from (3.10) we can see that structure 
of high spatial frequency appears naturally when the disordered state becomes 
unstable, and the critical wave vector (^-, ^) depends sensitively on system 
parameters. 



4. Dynamic Phase Transition and Pattern Formation 

By the dynamic transition theory [17, 18], the we know that as A crosses A c 
from below, the system always undergoes a dynamic transition to one of the three 
types, and the type of transition is dictated by the nonlinear interactions. In this 
section, we address the nonlinear interaction, the type and structure of the phase 
transitions, and the associated pattern formations. 

4.1. Transitions to LAM Patterns. In this section, we consider the case where 
the first eigenvalue is simple and the corresponding K\ which satisfies (3.6) is of 
the form K\ = (k\n/Li, 0, 0), k\ ^ 0. The type and structure of phase transitions 
of the system is dictated by the sign of the following parameter: 

As mentioned before, 

(4.2) lifxl 2 -^, A c ~2^. 

Hence the parameter B takes the following form: 

(4.3) 

V ; ' 3 27 y/a 

Theorem 4.1. When K\ = (kin/Li, 0, 0), k\ ^ Q, is the only wave vector which 
satisfies (3.6), the following assertions hold true: 

(1) If B < 0, then the phase transition of (2.13) at A c is Type-II, where 
\ c = (|i^i| 4 + cr)/|ifi| 2 and \K±\ 2 = klit 2 /L\. In particular, the prob- 
lem bifurcates on the side A < A c to two non- degenerate saddle points, and 
there are two saddle-node bifurcations at some A* < A c as shown in Figure 
1(a). 

(2) If B > 0, then the transition is Type-I, and the problem bifurcates on A > A c 
to two attractors U\ and u 2 as shown in Figure 1(b), which are given by 
(see Figure 2 for a structure of U\^): 

(4.4) u h2 = ±^^cob(*^) + o(\0 Kl \ 1/2 ), 

where is given in (3.5). 
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Figure 1: (a) Type-II transition given by Theorem ^.1: The disordered state u = 
is globally stable when A < A*, is metastable when A* < A < A c , and becomes 
unstable when A > A c . There are two saddle-node bifurcations at A*, (b) Tyep-I 
transition as given by Theorem 4-1- 



Figure 2: LAM pattern: a schematic structure of U1.2 given in (4.4). 



Two remarks are now in order. 

First, as we know, a Type-II transition correspond to the first-order transition 
in the sense of Ehrenfest. In this case, A* and A c in Figure 1(a) represent respec- 
tively the so called binodal and spinodal points. From (3.6), we know that A c is 
approximately 2y/cr, then by (2.9) and (2.12), we have the following formula for the 
spinodal which involves only \N and a as independent parameters: 

xN ~ ^ a 2 + (l-a) 2 



2{a{l-a)fl 2 2a 2 (l-a) 2 ' 

To our knowledge, there is no such simple formula for the binodal point. Nev- 
ertheless, some efforts have been devoted to figuring out A* numerically e.g. in 
[5]. 

In addition, when the control parameter A is between A* and A c , the disordered 
state u = is metastable. Perturbation may lead the system to other metastable 
states which may be far away from the original disordered one. 
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Second, Type-I transition corresponds to the second-order or higher-order tran- 
sition in physics. In this case the transition states are determined by the perturba- 
tion of the first eigenvectors. In Theorem 4.1, the transition states are Ui^, which 
are stable states, and represent a lamellar pattern as shown in Figure 2. 

4.2. Transitions to HPC Patterns. Now, we consider the case when the domain 
size satisfies the following condition: 

2 

(4.5) Li = 2ttL, L 2 = —j=TtL, and L 3 = 9irL. 

V3 

Here L and 9 are some positive constants, which are chosen in a way such that 
Kf := (f ,0,0) and K% := (^, ^,0) are the only wave vectors satisfying (3.6), 
where n is a positive integer. 

Theorem 4.2. Assume that the size of the domain satisfies (4-5) and K{ = 
(£, 0,0) and KZ, = (jj-, ^^,0) are the only two wave vectors which satisfy (3.6). 
Let B be the parameter defined in (4-1) with \K\\ — \Kf\ — \K%\. The following 
assertions hold true: 

(1) V 

72 = 0, 

then the phase transition of (2.13) at A c is Type-I, where A c is given by 
(3.6). The problem bifurcates on the side A > A c to an attractor which 
is homeomorphic to the one- dimensional unit sphere S . Y,\ contains eight 
non- degenerate steady states, with four saddle points W\, w 2 , w 3 andw^ and 
four minimal attractors U\,U2, U3 and W4 as shown in Figure 3(a). More- 
over, we have the following approximation formulas for the steady states: 

Ui l8 =±^P COB^) + 0(|A| 1/a ), 

-2,4 = ± J-^ 008(^*0 cos(^x 2 ) + am 1 ' 2 ), 

(4-6) 

Wl, 2,3,4 = ± 

where b and d are given in (5.16) and /?i(A) = Pk" is as in (3.5). 

(2) // 

72 ^ and B < 0, 

then the problem bifurcates on both sides of X c and the transition is Type-II. 
Moreover, there are four steady states bifurcated out on the side A < A C; 
including three saddle points and one unstable node. On the side A > A C; 
the problem bifurcates to two steady states, which are saddles. 

(3) // 

72 ^ and B > 0, 

then the transition is Type-Ill. Again, there are bifurcations on both sides 
of A c = 2-y/cr. On the side A < X c , there are two saddles bifurcating out 
from the origin. On the side A > \ c , the problem bifurcates to four steady 
states, including one stable node and three saddles. Moreover, there is a 



I /81(A) ,n 

y-bTTc coa{ d xi) 
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Figure 3: (a) For 72 = 0, the attractor after the transition is homcomorphic 
to S 1 and ui, 1*2, «3 and U4 are minimal attractors. (b) For 72 > 0, the sectorial 
regions V\ and V2 are separated by the lines 0u>i and 0t«4, and u\ is the minimal 
attractor in Vi. (c) For 72 < 0, the sectorial regions V\ and V2 are separated by 
the lines 0u>2 and OW3, and U3 is the minimal attractor in region V^. 



neighborhood V C H of u = 0, which can be decomposed into two disjoint 
sectorial regions V\ and V2 such that V — V\ U V 2 and the phase transition 
is first order in V\ and is n-th order in V2 with n > 2 in the Ehrenfest 
sense. V\ and V2 satisfy the relation given in (5.18). In region V2, there 
is exactly one minimal attractor as is shown in Figure 3(b) for 72 > and 
in Figure 3(c) for 72 < 0. The minimal attractors u\ and W3 can again be 
approximated as in (4-6). 

We note that, among all the steady states given in (4.6), are related to LAM 
as shown in Figure 2, 101,2,3,4 lead to HPC given in Figure 4, and U2,a correspond 
to rectangular patterns as in Figure 5. In the case of Type-I transition, as the 
disordered state loses its stability, the transition happens in two stages. There 
is first a fast transition from the disordered state to the bifurcated attractor E>, 
then followed by a slow transition within the bifurcated structure. Depending on 
the initial perturbations, the patterns finally emerging from the transition may be 
cither lamellae or rectangles. 

In Type-Ill transition case, transition may happen on both sides of the critical 
point A c . On the side A < A c , the disordered state is metastable, and perturbations 
may lead the system to some other metastable states far away from the disordered 
state as discussion in last section. On the side A > A c , the disordered state is 
unstable, and the transition may be either first order or second order depending on 
the initial perturbations. When the perturbation leads to second order transition, 
there is first a fast transition to the local attractor composed with two steady states 
with the HPC structure, one steady state with LAM structure, and two heteroclinic 
orbits connecting them as shown in Figures 3(b) and 3(c). The pattern eventually 
goes to LAM. 

4.3. Transitions to Squares and Spheres. In this section, we return to the 
situation when the first eigenvalue is simple, and we study the case when the cor- 
responding Kx satisfying (3.6) is of the form Ki — (kin/Li, k2TT/L 2 , 0), fej 7^ 0, 
1 < i < 2, or K\ = (kin/Li, fc 2 7r/L2, k 3 ir/L 3 ), fc, 7^ 0, 1 < i < 3. For simplicity, 
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(a) 



® ® • ® I 

5 ® ® ® * 
® ® « ® | 

') ® ® * ® 

6 Ok a 

(b) 



Figure 4: (a) HPC structures determined by 1^1,2,3,4 m (4-6). (b) Top view of the 
structures given in part (a). 



we will give results for the case k\jL\ = k^jLi = k^/L^, but the general situation 
can be dealt with in the same way. With this assumption, the first critical eigenvec- 
tors will have square pattern when K\ — (kin/Li, fc 2 ' r /-t'2j 0), and a cubic pattern 
when K\ = (kiir/Li, ^2"" /L2, k^ir/ X3). We define the following parameters: 

= 16/ 1 \Ktf \ a 

2 73 9 \2\Ki\ 4 -a 2U\K 1 \ i -2\\K 1 \ 2 ~a) 12 ' 

(4.7) B 3 = 7 3 LJJ ' 



2\K^ + 6A|Xi| 2 - 9cr 46|ifi| 4 - 6A|is:i| 2 - 9a 



1 



2 



Using (4.2), we have 
(4.8) 



+ 0^01^14 i08A|i^i| 2 - 18cr ' 72 



o 152 7 2 2 



Theorem 4.3. ^4ssMme that K\ = (kiir/Li, fc 2 7r/.L2, 0), k\jL\ ~ ^2/^2 7^ 0, is 
the only wave vector which satisfies (3.6). Then the following assertions hold true: 

(1) If 

B 2 < 0, 

then the phase transition of (2.13) at A c is Type-II, where the critical value 
A c is given by (3.6). In particular, the problem bifurcates from (u, A) 
Mi. A c ) on the side A < A c to two non-degenerate saddle points, and there 
are two saddle-node bifurcations at some A* < A c . 

(2) // 

B 2 >0, 
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Figure 5: Rectangular patterns: 
a schematic structure of 
given in (4.9). 




Figure 6: Sphere patterns 
in face-centered-cubic lattices: 
a schematic structure of ui,2 
given in (4.10). 



the transition is Type-I, and the problem bifurcates on A > A c to two attrac- 
tors U\ and U2 , which can be expressed as ( see Figure 5 for the structure of 
u%, 2 ) 

, 16@Ki fhnx! ( k 2 irx 2 n , 1/2 , 

(4.9) Ul,2 = ±4/ r cos(— — ) cos(— — ) + o(\f3 Kl I 1 ), 

where /3k 1 is as in (3.5). 

Theorem 4.4. Assume that the only wave vector which satisfies (3.6) is K\ = 
(fci7r/Li, k 2 7r/L 2 , k^Tr/L^), ki/L\ = k 2 /L 2 = k^/L^ ^ 0. Then the following 
assertions hold true: 

(1) // 

B 3 <0, 

then the phase transition of (2.13) at \ c is Type-II, where the critical value 
A c is given by (3.6). In particular, the problem bifurcates on the side A < A c 
to two non- degenerate saddle points, and there are two saddle-node bifur- 
cations. 

(2) // 

B 3 >0, 

the transition is Type-I, and the problem bifurcates on A > A c to two attrac- 
tors u\ and u 2 , which can be expressed as ( see Figure 6 for the structure of 

u l,2) 

(a m> i / 64/3^ Mkx x ,k 2 nx 2 ,k 3 Trx 3 1/2 

(4.10) u h2 = ±J 27 ^ Bs cos{— —) cos(— — ) cos(— — ) + oQPki I 1 ), 

where /3k 1 is as in (3.5). 

5. Proofs of Main Theorems 

Proof of Theorem J^.l. The proof relies on the following center manifold reduction, 
which reduce equation (2.13) to a one-dimensional ODE. 

Let u = v + $(y, A), where v = j/e^i, and $ is the center manifold function which 
we will calculate later. Thanks to the classical center manifold theorem [12, 28, 18], 
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we have (v, = and $(y, A) = o(|y|). Now multiplying both sides of (2.13) by 
e Kl and integrating over 0, we obtain the reduced equation to the center manifold: 

(5-1) tt =MX)y -' 1 M\ g{y)+0W% 
where is the volume of f2, and 

(5.2) g(y) = G 2 (y) + G 3 (y) + G 23 (y), 

(5.3) G 2 (y) = 72 / v 2 e Kl dx, 

Jn 

(5.4) G 3 (y) = 73 / v 3 e Kl dx, 

Jn 

(5.5) G 23 (y) = 2 72 f v*(y, X)e Kl dx. 

Jn 

By direct computation, we know 

(5.6) G 2 = f n (ycos(kiirxi)/Li) 2 cos(ki7rxi/Li)dx = 0, 

(5.7) G 3 = J (ycos(fc 1 7ra; 1 /L 1 )) 3 cos(fc 1 7rxi/L 1 )da; = ^|r% 3 . 

To evaluate G 23 , we need to compute the center manifold function <&(y, A). We 
will use the following second order approximation formula of A) (see e.g. [17] 
for details): 

(5.8) -£ A <%, A) = P 2 G 2 (y, A) + 0(|/3(A)||y| 2 ) + o(\y\ 2 ). 

Here P 2 : H — > ii^ is the canonical projection with /f^ 1 being the subspace of i7 
spanned by all stable eigenfunctions of L\, C\ is the restriction of L\ to H 2 , and 
/3(A) is the first eigenvalue of L\. 

The above formula (5.8) implies that 

$(y,A)= ^ \)e K + o(\y\ 2 ), where 

Jl f A ..2, 

(2p Kl (\)-p K ){e K , e K )i 



$ K (y,\)= — / Av 2 e K da; 

"7/r Jn 



I2PK / 2 , 



{2p Kl (\)-/3 K )(e K , e K ) H Jn 
Note that 

/ e 2 Ki e K dx= f 1 + C ° s{2kl7TXl/Ll) e K dx = unless K = (2fc l7 r/L 1 , 0, 0). 

For if = (2fci7r/ii, 0, 0), we have: 

(5.9) $ K (y,A) = -2£W, 

^2^2/^2 

Where C = 14fcXAf-2Afc^V^- ( T- 

Thus we obtain the center manifold function as follows: 

(5.10) d>(y, A) = -2£ l2 y 2 cos(^^) + o(\y\ 2 ). 
Plugging (5.10) into (5.5), we find 

(5.11) G 23 = ~miiy 3 + o(\y\ 3 ). 
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Finally, plugging (5.2), (5.6), (5.7) and (5.11) into (5.1), wc can derive the fol- 
lowing reduced equation of (2.13): 

(5-12) f = fa (A), - ^ (^73 - ^J^JS + o(M 3 )- 

It is known from [17, 18] that the transition type of (2.13) at the critical point A c 
is completely determined by (5.12). Thus all the assertions of the theorem follow 
from (5.12) except the claim for the saddle-node bifurcation in Assertion (1), which 
we prove as follows. 

8\K I 2 

When 73 < , ' ,4 — — 7I, we know from (5.12) that there are two branches 
00 yc 

of saddle points bifurcating from (u, A) = (0, A c ). The result will follow if we can 
show that there exists A < A c such that the following two conditions hold: 

(i) u = is globally stable for any A < A. 

(ii) The two branches of saddle points are uniformly bounded when A < A < A c . 

To verify (i), we do an energy estimate. For any A < A = — 72/273, and u ^ in 
H, multiplying the right hand side of (2.13) by u, and integrating over f2, we get: 

/ -\Au\ 2 + A|Vu| 2 - 2 72 w|Vu| 2 - 37 3 u 2 |Vu| 2 - au 2 dx 
Jn 

< / |Vu| 2 (A + 2|7 2 u| - 373U 2 ) - au 2 dx 
Jn 

< [ \Vu\ 2 (\- l3 u 2 -2 l3 (\u\-^) 2 + ^-)-au 2 dx 
Jn 273 273 

< / |Vu| 2 (A+ -^-) -ou 2 dx 
Jn 2 73 

< 0. 



Therefore, when A < A, (2.13) has no non-trivial singular points in H. It is known 
that for any A £ M 1 , (2.13) possesses a global attractor. Thus condition (i) follows. 
Condition (ii) is a direct consequence of the existence of global attractors for any 
A e R 1 . The claim for saddle-node bifurcations is thus proven, and the proof is 
complete. □ 



Proof of Theorem 4-2. The eigenvectors of L\ corresponding to K\ and are 

ei = cos(-xi), 

<5 - 13) ,» , > 



Let u = v + &(v, A), where v = y\e\ +1/262 and $ is the center manifold function. 
Again, $ can be approximated to the second order via the formula given in (5.8). 
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In this case, we have 



(5.14) 



3l2\KZ\ 2 yi y 2 ,3n ) 

U\KW-2\\KW-v C ° S( C0S ^ 2 ) 



cos(— — x 2 ) + o(\y\ ), 



2&\Kf\ 4 -4\\K?\ 2 -4cr v d 

where |A"f| = |^| = n/L. 

The reduced equation to the center manifold is now given by 

% = Pi(X)yi l -^p^yl + by\ + cy,y\ + o(\ y \% 
(5.15) *t 4 

-jjT = ftWlfc - l^i C | 2 7 2 yiy 2 + *y 2 3 + ej/ 2 y 2 + o(|y| 3 ), 
where ft (A) = fo.(A) = foe (A) and 

, 2|^i 4 72 2 

U\K^ - 2X\K^\ 2 - a 4 73 ' 



3|^| 4 72 2 _ 3|^| 2 

' 14|JCf | 4 - 2X\K1\ 2 - 2a 4 73 ' 

3\Ktf>y* 9|i^| 2 



73, 



(5.16) 



28|/q c | 4 -4A| J Fq| 2 -2cr 28|iff| 4 - 4A|-RTf | 2 -4a 16 
= 3|^| 4 7 2 _ 3\Kf\\ 

' 7\Kf\* - - o 2 73 ' 

As before, we have l-ft^l 2 — \fa. Then when A is close to A c , the above coefficients 
can be written as 

6^ 72 2 -^73 + o(|A-A c |), 
C ^ 7 2 -^73 + o(|A-A c |), 

e ^^7 2 2 -^73 + o(|A-A c |). 

It is known that the transition of (5.15) near A = A c can be understood by studying 
the behavior of the equation when A is at the critical point A c = (\Kf | 4 + a)/\Kf | 2 . 
At A = A c , /?i(A) = 0, and (5.15) reads (up to 3rd order terms): 

dyi l^i| 2 7 2 2,13, 2 

-rr = - A V2 + b Vi + c yiV2i 

(5.17) dt 4 

^ = -|^ c | 2 7 2 2/i2/ 2 + ^ 2 3 + e y 2 2; 2 . 
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(a) (b) 



Figure 7: The topological structure of flows of (5.15), (a) for 72 = and /?i(A) = 0, 
and (b) for 72 = and /3 X (A) > 0. 



Note that by (5.16) we have b + 4c = Ad + e, from which we see that on the 
straight lines 1/2 = ±2t/i, the equation (5.17) satisfies that 

C ^ = ±2, for all y x such that - \K\?^y\ + (b + 4c) j/? ^ 0. 
ayi 

Hence, the straight lines y 2 = ±2yi are orbits of (5.17). Obviously, 1/2 = also 
consists of orbits of (5.17). If 72 7^ 0, all straight line orbits of (5.17) lie on one of 
these three lines. If 72 = 0, y\ = also consists of orbits of (5.17), and in this case 
all straight line orbits lie on one of four straight lines y% = ±2yi, y 2 — or yi = 0. 

For the case 72 = 0, the numbers b, c, d, and e are all less than zero by (2.11) 
and (5.16). It is easy to see that all the straight line orbits of (5.17) tend to the 
origin as shown in Figure 7(a), which implies that the origin (2/1,1/2) = (0,0) is 
locally asymptotically stable. Therefore, the transition of (5.15) is TypeT [18]. 
Omitting o(|y| 3 ) terms in (5.15), the bifurcated structure are shown in Figure 7(b). 
The origin becomes unstable, and there are eight other steady states emerging 
from the transition, among which four of them are stable and the rest are saddle 
points. The stability of the steady states can be analyzed by linearizing (5.15) 
about the corresponding steady state, and the results are shown in Figure 7(b). 
By the center manifold reduction, the bifurcated structure for the original equation 
(2.13) is a perturbation of the structure shown in Figure 7(b), and Assertion (1) is 
thus proved. 

When 72 ^ 0, the straight line orbits of (5.17) consist of y 2 = ±2yi and y 2 = 
as shown before. The phase diagram of (5.17) near the origin can now be easily 
determined and are shown in Figures 8(a) - 8(d). 

For 72 7^ and b > 0, a neighborhood of the origin is divided into six regions by 
the straight line orbits, and four of them are hyperbolic and the rest two regions 
are parabolic as shown in Figures 8(a) and 8(b). Since orbits starting from any of 
these regions are eventually repelled away from the origin, the transition is Type-II 
[17]. The number and type of equilibrium points bifurcated out for each case can 
be analyzed in the same fashion as in the proof for Assertion (1), thus Assertion 
(2) is proved. 

For 72 ^ and b < 0, a neighborhood of the origin is also divided into six regions 
by the straight line orbits, and four of them are hyperbolic and the rest two regions 
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(c) (d) 



Figure 8: The topological structure of flows of (5.17), (a) for 72 > and b > 0, (b) 
for 72 < and b > 0, (c) for 72 > and b < 0, and (b) for 72 < and b < 0. 




(a) (b) 



Figure 9: The topological structure of flows of (5.15) for A > A c , b < and (a) for 
72 > 0, (b) for 72 < 0. The two sectorial regions 0\ and O2 are seperated by the 
rays 0^2 and OZ4. 



are parabolic as shown in Figures 8(c) and 8(d). Here, orbits in the parabolic 
regions tend to the origin, and the transition is Type-Ill [17]. The flow structure 
of the reduced equation (5.15) when A > A c can also be obtained easily and the 
result after dropping o(|y| 3 ) terms is shown in Figure 9. From here, we see that 
there is a neighborhood 77 of the origin of the yij/2 plane, which can be decomposed 
into two disjoint regions 77i and 77% such that 77 = 77 1 U77 2 and in region A/2, 
there is exactly one stable node bifurcating out on the side A > A c . Actually, 77i 
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and M 2 are slight perturbations of 0\ and 2 respectively as shown in Figure 9. 
Corresponding to Af, there is a neigborhood V of the origin in space H, which can 
be decomposed into two disjoint regions V\ and V 2 such that V = V\ U V 2 . Here V\ 
and V 2 are chosen such that 

(5.18) {1 - P 2 )Vi = Afi, 1 < t < 2, 

where P 2 is the canonical projection as in (5.8). Thus, assertion (3) is proved and 
the proof is complete. □ 

The proofs of Theorem 4.3 and Theorem 4.4 are similar to that of Theorem 4.1, 
and the details are omitted here. But for the convenience of the readers, we present 
the center manifold functions and the reduced equations to the center manifold. 

Proof of Theorem 4-3. By similar calculation, we have the following center manifold 
function: 

A) = -£1722/ cos( — ) -&722/ cos( — ) 

Li L 2 

2 2fci7rxi 2k 2 irx 2 2 
-6722/ cos( — )cos( — ) + o(\y\ ). 

L\ L 2 

Here £1, £2 and £3 are as follows: 
6(A) 

6(A) 



14|^ n | 4 - 4\K U \ 2 \K 12 \ 2 - 2\K 12 \* + 2X(\K 12 \ 2 - \K n \ 2 ) -a' 

\K12l 

14|X 12 | 4 - 4|^ 11 | 2 |^ 12 | 2 - 2|X n | 4 + 2A(|X n | 2 - \K 12 \ 2 ) -a' 



Sdv ' 14|X 1 | 4 -2A|K 1 | 2 -a' 

£■2^2 £-2_-2 h 2 -K 2 h 2 -K 2 

where l^nl 2 ^^, \K 12 \ 2 = |^| 2 = + 



^1 ^2 ^1 ^2 

The reduced equations of (2.13) to the center manifold in this case is: 



dy 
dt 



K1 WV- l^i| 2 (^73 - (6(A) + 6(A) + ^ 3 (A)) 72 2 )y 3 + o(|y| 3 ). 



Proof of Theorem 4-4- The center manifold function for this case is: 

«(2/, A) = -r/1722/ cos( — ) -772722/ cos( — ) 

^1 ^2 

2 ,2fc 3 7ra;3, 2 ,2k\-KX\ x ,2k 2 nx 2 . 
- 773722/ cos( — ) - 774722/ cos( — ) cos( — ) 

iv3 L\ L 2 

2 ,2fci7ra;i. .2fc 3 7rx3. 2 / 2fc 2 7ra;2, ,2fc 3 7ra;3, 

-%722/ COS( )C0S( )- 77672J/ cos( — )cos( — 1 

L\ L 3 L 2 L-j, 

2 2^17X1 2k 2 T7X 2 2k 3 7TX 3 2 

-777722/ cos( — )cos( — )cos( — )+o(\y\ ), 

Li\ Li 2 1^3 

and the reduced equations of (2.13) to the center manifold is given by: 

f t = P Kl (X)y l^il 2 (^73 - vil)y 3 + o(\y\ 3 ), 



□ 
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Here r\ and r\i (1 < i < 7) are given by 

fc?7T 2 /£? 



where 



^i(A) 
%(A) 

»?4(A) 

%(A) 
%(A) 
^(A) 



2[-2|^| 4 + 16 |^ + \{2\K X \ 2 - - a] 



2[-2\K 1 \ 4 + 16 ^- + X(2\K^ 4 |^) - a] 

2^2 / r2 



2[-2\ Kl \* + 16 |^ + A(2|^ - - a] ' 

\K12l 

2[-2\K 1 \ 4 + 16\K 12 \ 4 + 2\(\K 1 \ 2 - 2\K 12 \ 2 ) - a] ' 

\K13f 

2[-2\K 1 \ 4 + 16\K 13 \ 4 + 2X(\K 1 \ 2 - 2\K 13 \ 2 ) - a] ' 

\K23f 

2[-2\K 1 \ 4 + 16\K 23 \ 4 + 2\(\K 1 \ 2 - 2\K 23 \ 2 ) - a] ' 

l^il 2 

U\K X \ 4 -2\\K l \ 2 -a 



vW = vi + V2 + V3 + ^(vi + m + m) + 



L.2 2 pj 7,2 2 7,2 2 



3 

7,2 2 7,2 2 7.2 2 7.2 2 7.2 2 

1 2 _ K '2 n , ^3^ I jv- |2 _ ^'l 7 ' , ^2^ , ^S 71 

r2 + r2 ' I""- 1 ! ~~ r2 r2 + 7-2 



TZ 1 r2 ' I---M 72 ' 72 ' 7 

lv 2 -^3 ^1 ^2 ^ 



3 



□ 
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